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Small quantum fluctuations in solitons described by the cubic-quintic nonlinear Schrodinger equa- 
tion (CQNLSE) are studied with the linear approximation. The cases of both self-defocusing and 
self-focusing quintic term are considered (in the latter case, solitons may be effectively stable, de- 
spite the possibility of collapse). The numerically implemented back-propagation method is used 
to calculate the optimal squeezing ratio for the quantum fluctuations vs. the propagation distance. 
In the case of the self-defocusing quintic nonlinearity, opposite signs in front of the cubic and quin- 
tic terms make the fluctuations around bistable pairs of solitons (which have different energies for 
the same width) totally different. The fluctuations around nonstationary Gaussian pulses in the 
CQNLSE model are studied too. 
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I. INTRODUCTION 

It is well known that generalized nonlinear Schrodinger 
equation (NLSE) with a saturable response can sup- 
port multi-stable temporal or spatial solitons in one- and 
multi-dimensional cases Q, 0, 0, S El • The saturation 
of the Kerr nonlinearity may also prevent collapse of opti- 
cal spatio-temporal solitons ("light buttes"), in two- and 
three-dimensional cases 0, @ ■ The general form of the 
NLSE with saturable nonlinearity is (in the temporal do- 
main) 

iU z + U tt +T(\U\ 2 )U = (1) 

where U(z,t) is the local amplitude of the electromag- 
netic wave, z and t are, as usual, the propagation dis- 
tance and reduced time, and the function JF(|[/| 2 ) de- 
scribes the saturable nonlinear response of the medium. 
The simplest form of the latter function corresponds to 
the cubic-quintic (CQ) nonlinearity, with the cubic and 
quintic terms accounting for the self-focusing and self- 
defocusing, F(\U\ 2 ) = \U\ 2 — b\U\ A , where 6 is the ratio of 
cubic and quintic nonlinearity. Despite the obvious dan- 
ger of the spatio-temporal collapse, we will demonstrate 
that reasonable results may be obtained in the latter case 
too. 

The corresponding cubic-quintic NLSE (to be abbre- 
viated CQNLSE) describes light propagation in chalco- 
genide glasses and organic media , if nonlinear ab- 
sorption may be neglected. In fact, CQNLSE is the sim- 
plest model that makes it possible to study quantum fluc- 
tuations around bistable solitons (and non-soliton Gaus- 
sian pulses) in systems with saturable nonlinearity, there- 
fore results reported below may apply to a broader class 
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of optical media than those which are directly modeled 
by the combination of cubic and quintic terms (although 
stability of solitons is not the same in models with dif- 
ferent forms of the saturable nonlinearity [ll|). 

The theory of quantum fluctuations around optical 
solitons has been developed during the past 15 years by 
means of several techniques, which were applie d to the 
family of solitons in the NLSE model [j^lla . Il4| , higher- 
order solitons in the same system [l5j, self-induced- 
transparency solitons 0| , and Bragg solitons 0] . How- 
ever, quantum theory has not yet been worked out 
for bistable solitons, which we aim to carry out in 
this work, using the CQNLSE model in the simplest 
(l+l)-dimensional c ase, and applying the known back- 
propagation method [hS| . By means of this method, we 
investigate quantum fluctuations as a function of the dis- 
tance passed by the pulse. A measurement scheme which 
should enable the observation of the predicted effects is 
the homodyne detection, which measures the inner prod- 
uct between the given pulse and a local-oscillator gauge 
pulse, according to the projection interpretation of the 
homodyne detection |19|. 

First of all, we find that the optimal squeezing ratio 
of the quantum fluctuations around the pulse improves 
when both the cubic and quintic terms are of the focus- 
ing type; on the other hand, the squeezing ratio degrades 
when the cubic term is focusing and the quintic one is 
defocusing. The quantum fluctuations around two soli- 
tons with equal widths, belonging to a bistable pair, are 
found to be totally different due to the effects of the quin- 
tic nonlinearity. Quantum fluctuations around bistable 
Gaussian pulses (nonstationary ones) will also be studied 
to compare the results with those for the solitons. 

This paper is organized as follows: in section [HJ we 
derive the quantum CQNLSE model and the correspond- 
ing adjoint linearized equation for quantum fluctuations. 
The effect of quintic nonlinearity on quantum fluctua- 
tions, and the differences between equal-width solitons 
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FIG. 1: The amplitude-pulsewidth relation for the bistable 
solitons as given by Eq. 0}. Two pairs of bistable solitons 
with equal pulsewidth, A\.2 and Bi,2, are selected as examples 
for detailed investigation. 



belonging to a bistable pair, are highlighted in section 
lllll Section HVI deals with the fluctuations around non- 
stationary Gaussian pulses, and comparison with the case 
of the solitons. The paper is concluded by section Ivl 



II. QUANTUM CUBIC-QUINTIC NONLINEAR 
SCHRODINGER EQUATION 

The cubic-quintic version of the NLSE QJ is 

iU z + U tt + 2 x \U\ 2 U + 3-/\U\ 4 U = 0, (2) 

where the cubic coefficient is normalized to be x — ^l- 
In the case ofv > 0, exact solutions for bistable solitons 
are available 



U(z,t) = 



2/3 



J0z 



Vl + 4 7 /3cosh (2y/]h) + x 



(3) 



where (3 is a propagation constant (intrinsic param- 
eter of the soliton family), subject to the condition 
1 + 4 7 /3 > 0. The peak power of the solution 
A 2 = 2/3/ (VI + 47/3 + x) , is related to the full- width at 
half-maximum pulse width , r, by the following relation: 



cosh (qSt^JxA 2 + 7-4 4 



3x + 4 7 A 2 
2 7 ,4 2 



X - 



(4) 



The amplitude-pulsewidth relation for the bistable 
solitons, corresponding to Eq. Q, is displayed in Fig. 
^for x = +1 an( i 7 < 0. As is seen, the pulsewidth of 
the soliton cannot be smaller than a minimum (critical) 
value « 3.0 for 7 = -0.1 and « 3.42 for 7 = -0.13. 
Two pairs of bistable solitons with identical pulsewidths, 



marked as A (which is close to the turning point) and B 
(taken farther from the turning point), are chosen to be 
the illustrative examples. In particular, the pulsewidth 
of the bistable solitons belonging to Pair B is 3.5, and 
their amplitudes are 1.147 and 2.033, respectively. Soli- 
tons belonging to both branches of the curve in Fig. ^ 
relation are stable as classical solutions. 

In the quantum theory, the classical CQNLSE is re- 
placed by its quantum counterpart, with U(z,t) directly 
replaced by operator field variable U(z,t): 



U z = i(j tt + 2i X U^UU + ii^U^UUU. 



(5) 



The quantized field must satisfy the equal-coordinate 
Bosonic commutation relations, 



U(z,t 1 ),W(z,t 2 ) 
U(z,t 1 ),U(z,t 2 



U\z,t l ),U\z,t2) 



0(6) 



Equation JSJ, which provides for the Heiscnberg-picture 
description, can be derived from the Hamiltonian 



H 



so that Eq. 10 is tantamount to iU z — [U, H] (the 
replacement of the ordinary equal-time correlation rela- 
tions for quantum fields by the equal-coordinate ones can 
be justified, in the present physical context, within the 
framework of the canonical approach, as shown in Ref. 
|2S). 

Next, we substitute U = Uq + u into Eq. to 
linearize the equation around the classical solution Uq 
for the soliton containing a very large number of pho- 
tons. The linearized equation for the quantum fluctua- 
tions reads 



iutt + Aix\U a \ 2 u + 9i-f\U \' l u 
V u, 



(J) 



where V is an effective evolution operator. The operator 
perturbation field u satisfies the same equal-coordinate 
Bosonic commutation relation © as the unperturbed 
field operator U. 

To describe the quantum fluctuations, we need to find 
the corresponding adjoint field which satisfies the condi- 
tion 



(u A \Vu) - (V A u A \u), 
where the inner product is defined by 



(f\g) = 



4-oc 



(f*9 + /5 f ) dt. 



(8) 



(9) 



This definition of the inner product conforms to the prin- 
ciple that any physical observable can be expressed as 



3 



the inner product between a characteristic measurement 
function and the quantum-field operator of the perturba- 
tion The adjoint field which implements the condi- 
tion © obeys the following linear operator equation: 

u A = iuf t + Ai X \U \' z u A + 9ij\Uo\' i u A 
- 2i x U$u A * - 6ijU$Uoii A * 
= -V J 



~ >A u A . 



(10) 

It can be checked that the inner product © between so- 
lutions of Eqs. (7J and (|10|) is preserved in the evolution. 
Using this invariance, one can express the inner product, 
between the quantum perturbation field and a properly 
chosen projection function, at an output point, z = L, in 
terms of the input quantum perturbation at the initial 
point z — 0, 

(u A (z = L,t)\u(L,t)) = (u A (0,t)\u(0,t)), 

which is the basis of the back-propagation method. 

Then, it is easy to calculate the quantum uncertainty 
of the output field, knowing the statistics of the input 
quantum- field operators. In particular, the squeezing ra- 
tio of an observable / is calculated as 

vax[(f L (t)\u(L,t)] _ var[(/ o (t)|u(0,t)] 



R(L) 



(11) 



var[</z(t)|fi(0,i)] var[</z(t)|fi(0,t)]' 

where var[-] means the variance, /z,(t) is the projection 
function at the output point, and fo(t) is the back- 
propagated projection function. In the case of the ho- 
modyne detection, the output is the inner product of 
the input field operator with the local oscillator state 
|l9l l2lj . In the usual squeezing experiment, the mea- 
surement function fh{t) is a local oscillator pulse with 
following expression [recall Uo{z,t) is the classical solu- 
tion for the field pulse (in particular, soliton)], 



h(t) 



U (L,t)e 



it) 



y/f^dt\U (L,t)\* 



(12) 



where 8 is an adjustable phase shift between the local 
oscillator and the signal pulse of the homodyne detection. 
The optimal (minimum) vale of the squeezing ratio R{T) 
can be chosen, varying the parameter 9. 

Finally, we assume the input quantum perturbation- 
field operator corresponds to a coherent state, which is 
good assumption in most cases. Based on the formulation 
given above, in the next section we calculate the optimal 
squeezing ratios for solitons in the quantum CQNLSE 
model. 



III. QUANTUM FLUCTUATIONS ABOUT 
SOLITONS 



Effects of quintic nonlinearity on quantum 
fluctuations 



To start the analysis of quantum fluctuations around 
the soliton, we fix the self-focusing cubic nonlinearity co- 
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FIG. 2: The optimal squeezing ratio as a function of prop- 
agation distance for different values of the strength 7 of the 
quintic nonlinearity, while the cubic-nonlinearity coefficient is 
fixed, % — 1-0- When 7 > 0, squeezing is stronger, while for 
7 < 0, it is weaker. 



efficient to be x = +1 and vary the quintic nonlinear- 
ity coefficient 7, in order to study effects of the quintic 
nonlinearity on quantum fluctuations. The initial pulse 
is taken in the form of the soliton solution in equation 
@. The optimal squeezing ratio vs. the propagation 
distance is shown in Fig. For 7 = 0, the CQNLSE 
is reduced to the usual cubic NLSE, for which quantum 
fluctuations have been studied in detail [Ig. If 7 > 0, 
i.e., both the cubic and quintic nonlinearities are focus- 
ing, the optimal squeezing ratio is almost the same as 
the case of 7 = 0, provided that the propagation dis- 
tance is short enough. That is because the quintic non- 
linearity is not strong enough to affect CQNLSE soliton 
in a such short distance. However, for longer distances 
(> 5 soliton periods), the effect of the quintic nonlinear- 
ity accumulates, improving the optimal squeezing ratio. 
On the other hand, the quintic term with 7 < is de- 
focusing, causing degradation of the optimal squeezing 
ratio, due to partial compensation between the focusing 
(cubic) and defocusing (quintic) nonlinearities; if the de- 
focusing effect is strong (7 = 0.2), the optimal squeezing 
ratio quickly decreases even after propagating a short 
distance. 

Figure |3] shows the dependence of the squeezing ratio 
on the phase of the local oscillator in the homodyne de- 
tection when the quintic nonlinearity is defocusing 7 < 0. 
This dependence indicates quadrature squeezing of the 
quantum fluctuations around the CQNLSE solitons, like 
the case of the ordinary NLSE solitons (7 = 0). With the 
self- focusing quintic term, 7 > 0, the CQNLSE solitons 
also undergo quadrature squeezing. 
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FIG. 3: The squeezing ratio as a function of the phase of the 
local oscillator 9 used for the homodyne detection. In this 
case, x — +1.0, 7 < 0. 

B. Bistable solitons 

Here, we aim to compare the evolution of quantum 
fluctuations for bistable solitons belonging to the two 
pairs marked in Fig. Recall that the pair A is taken 
at 7 = —0.13, which is close to its turning point, 7 = 7 cr , 
and pair B is taken farther from its turning point, at 
7 = —0.1; the amplitudes of the solitons are, respectively, 
A 1 = 1.3045, A 2 = 1.5531, and B x = 1.147, B 2 = 2.033. 
The optimal squeezing ratios for these two pairs of the 
bistable solitons are shown in Fig. 

We observe that, for the solitons belonging to the pair 
A, which have the same widths, r = 3.5, and slightly 
different amplitudes, the evolution of the quantum fluc- 
tuation is similar. One may expect that the soliton with 
a higher amplitude will be more squeezed, this is why 
the optimal squeezing ratio for the A 2 soliton is smaller 
than for the Aione, in the beginning. But due to the de- 
focusing effect of the quintic nonlinearity (recall we now 
consider the case of 7 < 0), the optimal squeezing ra- 
tios for both solitons degrade after propagating a certain 
distance (about 3 and 4.5 soliton periods for A 2 and A\, 
respectively), which leads to crossing of the squeezing- 
ratio curves for these bistable solitons. 

For the pair B, the difference of the amplitudes in 
the soliton pair is large. In accordance with this, we 
find that properties of the quantum fluctuations are to- 
tally different for these two solitons. The curve of the 
squeezing ratio for the soliton with the smaller ampli- 
tude, B\ = 1.147, is like the ordinary NLSE solitons 
(the dashed line in Fig. 0J, while the squeezing for the 
soliton with the larger amplitude, B 2 — 2.033, is much 
poorer, which is naturally explained by the strong defo- 
cusing effect exerted on the latter soliton by the quintic 
self-defocusing term. 
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FIG. 4: The optimal squeezing ratios for two pairs of bistable 
solitons. The case of the cubic nonlinear Schrfidinger equa- 
tion, 7 = 0, with the amplitude A = 1.0 is also plotted 
(dashed line) for comparison. 



IV. QUANTUM FLUCTUATIONS ABOUT 
NONSTATIONARY PULSES 

For situations where exact soliton solutions for pulses 
are not available, the variational approximation (VA) is 
known to be an efficient analytical method (see a review 
|2^1. In particular, the VA for the classical CQNLSE 
solitons was developed in Ref. [2i|, using the Gaussian 
ansatz for the pulse waveform, 

U(z, t) = A exp (- (t 2 /2a 2 ) + iat 2 ) , 

where a and a are, respectively, the width and chirp of 
the pulse. The use of the Gaussian makes sense not only 
because it is convenient for the application of VA, but 
also due to the fact that laser sources usually produce 
pulses with this shape (including the intrinsic chirp). 
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FIG. 5: The optimal squeezing ratio vs. the propagation 
distance for different nonstationary Gaussian pulses with dif- 
ferent widths. The result for the NCQSE soliton is plotted 
by the dashed-dotted line, for comparison. In this figure, 
7 = +0.1, and A = 1.0. 

Strictly speaking, the Gaussian ansatz may only pro- 
duce a nonstationary pulse (plus some radiation). Nev- 
ertheless, it is possible to calculate quantum fluctuation 
around it, and compare the results with those presented 
above for the solitons. In general, the squeezing ratio for 
nonstationary pulses cannot be larger than for the soliton 
of the same width, due to emission of radiation waves by 
the nonstationary pulse. Results of the calculation of the 
optimal squeezing ratio for the Gaussian pulses with dif- 
ferent widths are displayed in Fig. |SJ The Gaussian pulse 
produces a similar but poorer squeezing curve than the 
soliton of the same width, a = 1.25 (the dashed-dotted 
line, marked by "7 = 0.1" in Fig. |SJ). Better squeezing 
ratios can be obtained on a short propagation distance 
for shorter pulses, with a < 1.25, but they degrade very 
quickly with the increase of the distance. For broader 
pulses, with a > 1.25, the radiation modes will strongly 
affect the optimal squeezing ratio from the very begin- 
ning. 

When 7 < 0, the VA predicts bistable Gaussian pulses 
for any energy E — a\A\ 2 [23[. Figure © shows the 
optimal squeezing ratios for a pair of the thus defined 
bistable nonstationary Gaussian pulses, and they are 
compared with the curve for the NLSE solitons (marked 
by "7 = 0" ) . This pair of the bistable Gaussian pulses 
have a common width, a — 1.3, but drastically different 
energies, Eq — 1.30 and Eq — 6.66, respectively. Like in 
the case of the bistable solitons belonging to the pair B 



in Fig. 01 the pulse with the smaller energy shows fluc- 
tuations similar to those of the NLSE soliton, while the 
pulse with the larger energy has a much poorer squeezing 
ratio. 
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FIG. 6: The optimal squeezing ratios for nonstationary 
bistable Gaussian pulses, in the case of 7 = —0.1 and a = 1.3. 
The squeezing ratio for the ordinary nonlinear-Schrodinger 
soliton (the dash-dotted line, 7 = 0.0) is plotted for compari- 
son. 



V. CONCLUSION 

In this work, we have applied the back-propagation 
method to study the quantum fluctuations around 
bistable solitons described by the cubic-quintic nonlinear 
Schrodinger equation. It was found that the squeezing 
ratio for the soliton strongly depends on the sign of the 
quintic nonlincarity: the self-focusing quintic term helps 
to squeeze the fluctuations, while the self-dcfocusing one 
makes the squeezing poorer. In particular, the quantum 
fluctuations around bistable solitons seem totally differ- 
ent due to the self-defocusing quintic nonlinearity, al- 
though the solitons have the same width. Quantum fluc- 
tuation about nonstationary Gaussian pulses were also 
explored; in that case, radiation loss makes the squeez- 
ing weaker. 
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